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It is shown that the endpoint of the first order transition line which merges into a crossover 
regime in the phase diagram of the Nambu-Jona-Lasinio model, extended to include the six-quark 
't Hooft and eight-quark interaction Lagrangians, is pushed towards vanishing chemical potential 
and higher temperatures with increasing strength of the OZI-violating eight-quark interactions. We 
clarify the connection between the location of the endpoint in the phase diagram and the mechanism 
of chiral symmetry breaking at the quark level. We show how the 8q interactions affect the number 
of effective quark degrees of freedom. We are able to obtain the correct asymptotics for this number 
at large temperatures by using the Pauli-Villars regularization. 
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I. INTRODUCTION 



The last two decades have witnessed great efforts towards the understanding of the QCD phase diagram, in terms 
of effective low energy theories paralleled by QCD lattice calculations, see e.g. the recent reviews [H-[I|]> or the paper 
The domain of small to moderate baryonic chemical potential < \i < 400 MeV and temperatures < T < 200 
MeV, is of specific relevance for relativistic heavy ion collisions. Of fundamental importance in the study of the phase 
diagram are chiral symmetry and confinement, however the finite size {m u ,md,m s ^ 0) as well as the difference in 
the bare quark masses (m u =/= rrid m s ) pose major problems both from the calculational point of view and in the 
implications due to deviations from the ideal situations, where the chiral condensate and the Polyakov loop are known 
to be the appropriate order parameters to characterize the phase state of the quark-gluon system. 

The present study focuses on the chiral symmetry breaking aspects related to non-zero current quark mass values. 
Our arguments will be based on the successful model of Nambu-Jona-Lasinio (NJL) 6], combined with the Ua(1) 
breaking 2Nf flavor determinant of 't Hooft Hi- HI (NJLH). Moreover the most general U(3)l <8> U(3)r chiral invariant 
non derivative eight-quark (8g) interactions [l(| are included. These terms were proven to render the effective potential 
of the NJLH model globally stable, and their effect has been thoroughly studied in low energy characteristics of 
pseudoscalar and scalar mesons , at finite temperature [IH, [l3[ and in presence of a constant magnetic field [ljj ■ 
These studies have lead at instances to sizeable and unforeseen effects. Of particular importance for the present work 
is that the strength of the 8q coupling is strongly correlated with the temperature and slope at which the crossover 
occurs and that it can be regulated together with the four-quark coupling, leaving the meson spectra at T = 
unaffected (with exception of the scalar a meson mass which decreases with increasing 8q coupling) [lj| . As a result 
the symmetry breaking for large 8q couplings is induced by the 6q 't Hooft coupling strength, as opposed to the case 
with small 8q coupling, where the dynamical breaking of symmetry is controlled by the 4g coupling strength [Tlj |. We 
would like to comment on a natural question which arises here, namely, can higher order many-quark interactions be 
also important? With regard to this, explicit argumets of A. A. Andrianov and V. A. Andrianov are known [15| . which 
show that the structure of the QCD-motivated models at low energies with effective multi-fcrmion interactions and a 
finite cut-off in the chiral symmetry-breaking regime should contain only the vertices with four, six and eight-fermion 
interactions in four dimensions. This result explains partly the approximation used in our work. 

Thus, in this paper we give a quantitative account of local multi-fermion forces on the phase diagram in the (T, /x) 
plane, by comparing the results for two sets of parameters, in the small and large coupling regimes of the 8q strengths, 
corresponding to the two above mentioned alternative mechanisms of chiral symmetry breaking. 

Throughout the paper we work for simplicity in the isospin limit m u — nid ^ m s , breaking explicitly the chiral 
SU(3)l ® SU(3)n symmetry to the SU{2)i ® U(l)y (isospin-hypercharge) subgroup, and take the same baryonic 
chemical potential /i for all quark species. Generalizations to take into account the nonzero isospin chemical potential 
can be implemented as for instace in [Ui . 
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II. EFFECTIVE LAGRANGIAN 

The explicit form of the multi-quark Lagrangian considered is presented in [ToL [ill ] 

£ cf f = g(z7^ - m)q + C 4q + C 6q + 4J + 4? ■ (!) 
Quark fields g have color (N c = 3) and flavor (Nf — 3) indices which are suppressed, \i = 0, 1, 2, 3. Here 

A 9 - | [(9A a q) 2 + {qij 5 \ a q) 2 ] , (2) 
£ 6 <? = «(det qP L q + det gPflg), (3) 
4? = 8.91 mPB.qm){q m PLqi)] 2 , (4) 

A? = 16 52 [(qiPRq m )(q m PLqj){qjPRqk){qkPLqi)] • (5) 



The matrices acting in flavor space, A a , a = 0, 1, . . . , 8, are normalized such that tr(A a Ab) = 25 a b', Aq = yj 1 1, and 

Afc, k = 1,2,..., 8 are the standard SU(3) Gell-Mann matrices; Pl,r = (1 T7s)/2 are chiral projectors and the 
determinant is over flavor indices. The large N c behaviour of the model is reflected in the dimensionful coupling 
constants, [G] = M~ 2 , [k] = M~ 5 , [gi] = [g 2 ] = M~ s , which count as G — l/N c , k ~ 1/Nc f , and gi,g 2 ~ 1/A 4 or 
less. As a result the NJL interaction @ dominates over £ eq at large N c , as one would expect, because Zweig's rule 

is exact at N c = oo. Let us note that the 8q-interaction breaks Zweig's rule as well. 

Since the coupling constants G, k, gi,g 2 are dimensionful, the model is not renormalizable. We use the cut-off A to 
render quark loops finite. The global chiral SU(3)l x SU(3)r symmetry of the Lagrangian (JTJ at m = is sponta- 
neously broken to the SU(3) group, showing the dynamical instability of the fully symmetric solutions of the theory. 
In addition, the current quark mass m, being a diagonal matrix in flavor space with elements diag(m u , mj, m s ), explic- 
itly breaks this symmetry down, retaining only the reduced SU(2)j x U(l)y symmetries of isospin and hypercharge 
conservation, if m u = ma ^ m s . 

The model has been bosonized in the framework of functional integrals in the stationary phase approximation 
leading to the following effective mesonic Lagrangian £bos at T = // = 

£ c ff -> £bos = Cgt + £ q i, £ st = h a a a + 2 fl ab a a cr b + T^itVa^b + C(field 3 ), 

W ql (a,<f>) = -ln| detD^s | = - / ^L^i-Mh) = I d 4 x s £ ql , 

b = i, 6i = -y, & 2 -^ + ^|A 8 y, 

Y = ij a (d a a + i^d a 4>) + o 2 + {M, a} + <j> 2 + ij 5 [a + M,(f\ (6) 

written in terms of the scalar, a = A a cr a , and pseudoscalar, <f> = \ a (j> a , nonet valued quantum fields. The result of 
the stationary phase integration at leading order, £ st , is shown here as a series in growing powers of a and <f>. The 
result of the remaining Gaussian integration over the quark fields is given by W q \. Here the Laplacian in euclidean 

space-time D [ E D E = M 2 — d 2 + Y is associated with the euclidean Dirac operator D% = i-f a d a — M — a — 175^ (the 
7 Q , a — 1,2,3,4 are antihermitian and obey {70,7/3} = — 2(5 Q ^); M. — diag(M„, M<j, M s ) is the constituent quark 
mass matrix (to explore the properties of the spontaneously broken theory, we define quantum fields a a , 4>a as having 
vanishing vacuum expectation values in the asymmetric phase). 

The expression for the one-quark-loop action W n \ has been obtained by using a modified inverse mass expansion of 
the heat kernel associated to the given Laplacian [l7[. The procedure takes into account the differences A us = M 2 —M 2 
in the nonstrange and strange constituent quark masses in a chiral invariant way at each order of the expansion, 6, 
being the generalized Seeley-DeWitt coefficients of the new series. This modification distinguishes our calculation 
from the one made in (la |. In fact we consider the series up to and including the order b 2 that corresponds to the 
first nontrivial step in the expansion of the induced effective hadron Lagrangian at long distances. At this stage 
meson fields obtain their kinetic terms, but are still considered to be elementary objects. The information about their 
quark-antiquark origin enters only through the coefficients such as the average 

h = I [HMD + MM 2 ) + MM 2 )] (7) 
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over the 1-loop euclidean momentum integrals Ji with i + 1 vertices (i = 0, 1, 



^ 2 ) = 16^ + 1)/^PA^ +M2)m . 



(8) 



For the explicit evaluation of Ji(M 2 ) we use the Pauli-Villars regularization method with two subtractions in the 
integrand. The procedure is fully defined by the insertion of the particular operator 



p A = 1 - 1 - A 



d 



dM 2 



exp A 



d 



DM 2 



(9) 



Here the covariant cut-off A is a free dimensionful parameter which characterizes the scale of the chiral symmetry 
breaking in the effective model considered. To the order of the heat kernel series truncated, only the integrals Jo, J\ 
are needed. These are quadratic and logarithmic divergent respectively with A — » oo, all other J, are finite. Note 
that the recurrence relation 



J+i(M 2 



d 



dM- 



■MM 2 



(10) 



is fulfilled. If Ji is known for one value of i, then the function may be computed for other values of i by successive 
applications of the relation. 

In C st the h a are determined via the stationary phase conditions. These conditions and the pattern of explicit 
symmetry breaking show that in general h a can have only three non-zero components at most with indices a = 0, 3, 8, 
i.e. h a X a = diag(/i«, hd, h s ), which can be found from a system of three independent equations 



f Gh u 

Gh d 



■ A„ + — h d h s 
16 

A d + —z Kh s 
16 



91 
4 
91 



Gh s + A s + — h u h d 
16 



J hd{h 2 u 

| hs{hl 



hi 



hi 



h\ 



-hi) 
hi) 

hl)- 



9 4 hi = o. 



(ii) 



(1 2) 

Here Af = Mf — nif, f = u,d,s. The matrix valued constants of higher order, like for instance /iv in £ st , are 
uniquely determined once the hf are known [lll.[l9|. The stability of the effective potential is guaranteed if the system 
(ITU has only one real solution. For that the couplings must fulfill the inequalities g\ > 0, g\ + 3g2 > 0, Gg\ > 
(n/16) 2 . 



TABLE I. Parameters of the model: m u = m d) m s (MeV), G (GeV -2 ), A (MeV), k (GeV -5 ), g u g 2 (GeV -8 ). We also show 
the corresponding values of constituent quark masses M u = Md and M s (MeV). 





m u 


m s 


M u 


M s 


A 


G 




5i 


~92 


a 


5.9 


186 


359 


554 


851 


10.92 


1001 


1000* 


47 


b 


5.9 


186 


359 


554 


851 


7.03 


1001 


8000* 


47 



TABLE II. The masses, weak decay constants of light pseudoscalars (in MeV), the singlet-octet mixing angle 8 P (in degrees), 
and the quark condensates (uu) = (dd), (ss) expressed as usual by positive combinations in MeV. 







m K 


m v 


m v > 


u 


Jk 


dp 


— (uu) 3 


— (ss) 3 


a 


138* 


494* 


477 


958* 


92* 


117* 


-14 


235 


187 


b 


138* 


494* 


477 


958* 


92* 


117* 


-14 


235 


187 



TABLE III. The masses of the scalar nonet (in MeV), and the corresponding singlet-octet mixing angle 9 S (in degrees). 





m a (980) 


(soo) 


m /o(600) 


TO /o(1370) 


9 S 


a 


980* 


1201 


691 


1368 


23 


b 


980* 


1201 


463 


1350 


19 
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In this paper we use the two parameter sets of Table I, which differ only in the choice of the Aq coupling G and 
the 8q strength gi . Set (b) is the same as in [l3| (there was a misprint in the value for the constituent strange quark 
mass, which we corrected). Tables II-III display the numerical fits at T = /i — (input is denoted by a *). The only 
difference in the observables of the two sets occurs in the singlet-octet flavor mixing channel of the scalars, mainly 
in the a- meson (i.e. /o(600)) mass. The model parameters are kept unchanged in the calculation of the T and fx 
dependent solutions of the gap equations (see next sections). 

It is worthwhile to stress that there is an essential difference between the two alternative ground states chosen 
here as the configurations on top of which the T ^ and /i 7^ effects are studied: Case (a) corresponds to the 
standard picture of the NJL hadronic vacuum. In this picture chiral symmetry is spontaneously broken at T = fi = 
when G > G cr n. Case (b) corresponds to a new alternative, related to the pattern where G < G cr a. In this case 
chiral symmetry can be broken only due to the six-quark interactions, when \k\ exceeds some critical value (the 8q- 
interactions could in principle also induce symmetry breaking, however the mass spectra are then not well reproduced). 
One can hardly distinguish between the two cases at T = it, = 0, the spectra of v and ++ low-lying mesons do 
not show much difference: the model parameters are the same, except for the correlated G and g\ values. The larger 
value of g± in the case (b) is a signal of the increasing role played by the eight-quark OZI-violating interactions, but 
this does not affect the value of the mixing angle 9 p , and only slightly diminishes 9 S . Such insensitivity follows from 
the observation that the stationary phase equations (|11[) and mass formulae of the light ^ and ++ states 
only depend on the couplings G and g\ through the linear combination £ = G + gi(h^ + h\ + /i 2 )/4, except for the 
00, 08 and 88 states inside the scalar nonet. However as soon as T or /1 are finite, the hf start to change due to their 
intrinsic T, /1 dependence, acquired through the coupling to the quark loop integrals in the gap equations, eqs. ([12")) 
below. The T, /z dependence of the combination £ above is steered by the strength g\ . This is the main reason why 
the 8g-interactions may strongly affect the thermodynamic observables, without changing the spectra at T = fj, = 0. 



III. THERMODYNAMIC POTENTIAL 



A. Case of vanishing temperarure and chemical potential 

Before addressing the thermodynamical potential it is instructive to briefly discuss the effective potential of the 
model at T = fj. = 0. Using standard techniques [l9j . we obtain from the gap-equations 

hf + ^M f J (Mj) = (12) 

the effective potential V(Mf) as a function of three independent variables Mf = {M u ,Md,M s }. If the parameters of 
the model are fixed in such a way that eqs. (|11[) have only one real solution, the effective potential is 



V(M f ) = -- h f &M f --^ Yl M f J (M 2 f )dM f 



f—u,d,s f=u,d,s 



= i a)+^A+y(^) +¥}] +_1 ^ j_ l(M/ 2 ), (13) 



8tt 2 



f=u,d,s 

where h 2 = h\ + h\ + h 2 , hj — h\ + h\ + , and we extend definition (J7J for index i = — 1 with 

J_i(M 2 ) = - I" J (M 2 )dM 2 = -I (m 2 J (M 2 ) + A 4 In (l + ^] ) . (14) 



Here Jo has the explicit form 



2 ~uv^ A2 



Jo(M 2 )=A 2 -Af 2 ln(l + ^) (I:, 1 



for the given choice of regulator. 

The first integral in (fTB"]) accounts for the leading order stationary phase contribution. The second integral describes 
the quark one-loop part. Since both integrands in p^|) are exact differentials, the line integrals depend only on the 
end points. The low limit of the integrals is adjusted so that V(0) = (to understand this, it is enough to notice that 
the power-series expansion of V{Mf) at small Mf starts from V(0); this term does not depend on Mf and, therefore, 
does not affect the physical content of the theory; so we simply subtract it, calculating the potential energy of the 
system with regard to the energy of the symmetric vacuum in the imaginary world of massless quarks). 
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B. Case of finite temperature and chemical potential 

The extension to finite T and [i of the bosonized Lagrangian © is effected through the quark loop integrals Ji 
(see eq.©). Due to the recurrence relation (fTU|) it is sufficient to get it just for one of them, J , by introducing the 
Matsubara frequencies, w n , and the chemical potential, /i, through the substitutions [20l | 

oc 

/ dp 0E -> 2ttT ^2 > PoE^>u n -ifi, uj n = nT(2n + 1). (16) 

n= — OO 

Inserting (fT6j) into Jo we obtain 



J (Af 2 ) -> J (M\T,fi) = 16tt 2 T J E 



(2tt) 3 r f-^ C„ - 2ifj,TrT(2n + 1) 



C„ = £ 2 +^ 2 T 2 (2n+l)-- M 2 , E p = ^JhP+p 2 . (17) 
The sum over n is evaluated to give 

El 1 1 1 / fJL + E p fJL — E p \ 

_ C n - i2yu n = Wf? (n + a)(n + b) = WE p I™ ^T" " tanh ) (18) 

where we use the abbreviations a = h + (-Ep — A*) anc > ^ —\~ j^t (^p Thus we have 

J (M 2 ,7» = 4 / dl^HpE^A— (1 - n, - n ¥ ) = J (A/ 2 ) - 4 / d|^||p £ | 2 p A ^— -i 
Jo &p Jo ^P 

where the quark, anti-quark occupation numbers are given by 

1 1 



(19) 



1 + e t 1 + e t 



(20) 



Notice that J (Af 2 ,0,0) = Jq(M 2 ), therefore the vacuum piece is well isolated from the matter part. The remaining 
integral containing the quark number occupation densities n q ,riq is strictly finite, the A dependent terms being a 
remnant of the Pauli-Villars regularization scheme. At small T and M ^ we have 



J™d\p E \\p E \ 2 ^ „ J^\ M \M_ e -¥ = T V2TMe-^ J~*xy/5e-y[\ + g: 
Vhl 3 m ( 3T \ 

T^e^T i + + ... . (21 

2 V 8M J y ' 

The special feature of this integral is that it vanishes exponentially with T — » 0. 

Now we are ready to evaluate the thermodynamical potential. Indeed, the gap-equations at finite T and /i are 

h f + ^M f J o (M},T^) = 0. (22) 

Consequently, 

V(M f ,T, A*) = - 5 / E M M / " 74 / E M /^o(M 2 , T, M )dM/ + C(T, M ) 



f—u,d,s f—u,d,s 



-•■ f=u,d,s 



— UGh 2 + nh u h d h s + JL. (tf f ) + 3g 2 h f 



^ J2 J-i(M 2 f ,T,f,) + C(T,^), (23) 



where the function C(T, fi) does not depend on M , and therefore cannot be determined from the gap equation; it will 
be found from other arguments in the end of this section, but obviously C(0, 0) = 0. The integral J_i(Af 2 ,T, fi) is 
the immediate generalization of the T = /j, = case (|14p 



M 



J_i(M^, T,fx) = — / J (M , T, M )dM 2 = J_i(M 2 ) + J^M 2 , T, //), (24) 
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where the medium contribution to J_i(M 2 ,T, fx) is 

J™t\M\T,n)=l dM 2 d|p £ ||p i5 | 2 p A ^±^=4/ d\p E \\PE?PA PE dM 2 ^±^ 
Jo Jo Jo JO 



= 8^ d\p E \\p E \ 2 p A p E (2 (E P (M) - E p (0)) + Tin 

; )(> 



n q Mn-qM \ 
n q0 n q0 



f \ 

8T d\p E \\pE\ 2 PAp E l*-r 
Jo 1 



gp(Q)-M \ / EpjO) + M 

1 + e t 1 + e r 



(25) 



(l + e t J (l + e — — — j 

with n q o,n-qo and n q M,n q M referring to the occupation numbers for massless and massive particles correspondingly, 
E p (M) = y/M 2 +p E , E p (0) = \p E \. In spite of the fact that the integral J"\ cd (M 2 , T, (i) is convergent, we still 
keep the regularization p\ to be consistent. Note that the action of the operator /5a (see eq. ©) on any smooth 
function, depending on M 2 through the energy, f(E p (M)), can also be expressed in terms of momentum as pAf(E p ) = 
PA PE f(E p ), where 

Pa Pe = 1 - ( 1 - A 2 4l) e*P f A 2 i) • (26) 



dpi 



B 



Then, noting that, for instance, 



&We\\pe\ 2 Pa Pe In 



e 



|p B | 3 / gajMW. 
PA -^r~ In 1 + e 



3 





- [°° a\x 1 \PE\ 4 n qM 1 f°° A ^ ||- ,4- n <?Af 
PA i d|PE| 3T^(M) = 3Ti diP£i|P£i ^E^M 



^/ d| ^^-^| ln ( 1 + e 



where we used the fact the surface term disappears, we get finally 



(27) 



T med/u2 t \ 8 / A\" \\X |4~ j n qM + n qM n q0 + n qQ \ 

J_ x (M ,T )M ) = -- ^ d|p £ |b £ | p Afe [-^^ ^-J. (28) 

It is important to realize that the expansion of the integral (|28|) for small values of T and at = starts from the 
term 

ff(M 2 ,r,0) = ^ r lfe|3d| f g| + C(Tf Mle-T) = ±Vt 4 + 0(Tf Mle-"). (29) 
3 Jo e -¥- 45 

This leading contribution in T arises from the combination ~ (n q o +n q o) and does not depend on the cut-off, i.e. the 
Pauli-Villars regulator pA PE does not affect the leading order of the low-temperature asymptotics of the integral. To 
make this clear let us consider the typical integral in eq. 



<1\Pe\\pe\ = tttz = / e-titt^T / dx 7— TTm ■ (30) 



E P {M) J M 1 + e ^ Jo l + e*+f 

If M — 0, then the integral can be evaluated explicitly and is found to be -jfgT 4 , leading us to the result (|29|) . If 
M 7^ 0, we obtain at once the estimate at small T 



1 + e x +^ 4 \ T J 

We conclude that the integral vanishes exponentially for small T and does not contribute to the leading order term 
in eq. (f29|) . It is then clear that the action of the Pauli-Villars regulator pA PE on the integrand, which consists in 
subtracting the contribution of the massive Pauli-Villars states, will not affect the leading term as well. 

One might worry that the low-temperature expansion of the integral (|29p starts from the unphysical contribution 
which corresponds to the massless quark states. This fear would be valid if the potential had not the term C(T, //). 
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Let us fix the M/-independent function C(T, (j,) in eq. (|23|) to avoid the problem. For that it is instructive to compare 
the matter quark-loop part of the thermodynamic potential obtained here (i.e. the J™° d -part) with the corresponding 
result of the standard NJL approach. There is only one difference between such calculations: we use the Pauli-Villars 
subtractions instead of a 3-dimensional cut-off A3 . Therefore we can expect that if one removes the regularizations 
in both approaches {pAp E ~> 1 and A3 — > 00) the finite matter part of the thermodynamic potentials must coincide. 
From this requirement of consistency we find 

C(T ,„, . ^J~ dM «-^ a«r,o) . -1Ml&. ( 32 ) 

As a result, the unwanted massless quark contribution to the vacuum energy at small T disappears. 



C. Effective number of quark degrees of freedom 



A quantity of interest related with the thermodynamic potential is the number of quark degrees of freedom present 
at a certain temperature. For that we consider the quark pressure difference from the zero-temperature value 

_ p(T)- P (0) 

V{1) - ^ 2 T 4 /90 ' W 

in order that the total ^(0) = 0. Here p(T) = —V(M'f,T,fi — 0), where MJ denotes the gap equation solution at a 
given T . Dividing the pressure in eq. (I33[) by 7r 2 T 4 /90, the result is presented in Stefan-Boltzmann units, i.e. one 
explicitly counts the number of relevant degrees of freedom. 

It is known that in the case of massless quarks of two and three flavors the fermionic degrees of freedom at T > T c 
are estimated as v = (7/8) x3x2x4 = 21 and v — (7/8) x3x3x4 = 31.5, respectively. The system studied here 
consists of three-types of light quarks: u, d, s. Hence v is expected to be in the interval 21 < v < 31.5 in the region 
T > T c , where chiral symmetry is "restored" (up to the explicit symmetry breaking effects caused by the current 
quark masses). Indeed, the solid curves v(T/T c ), plotted for the sets (a) and (b) in fig. [TJ at some values T/T c > 1 
enter the interval and approach asymptotically the upper bound v — 31.5 at high T: in particular, we have already 
at T/T c = 2.5 that i/(2.5) = 30.95 (set (a)) and v(2.5) = 30.25 (set (b)). This is too fast as compared with the lattice 
estimates [2l| in 2+1 flavor QCD, but the difference can be ascribed to the simplifications introduced by the model 
under consideration (the essential difference is that the NJL model does not possess the quark-confinement property 
of QCD). 

We can gain some understanding of the asymptotic behavior of v{T) by considering eqs. (fT9| and (|24|) . Firstly, it 
is easy to see that the integral Jo(M 2 , T,^)->0atT^oo. The reason for this is very simple and is contained in the 
integrand (1 — n q — nq) which vanishes at T — > 00, as it explicitly follows from eqs. (|20[) . Secondly, eq. (|24l) contains 
Jo as an integrand. Thus we conclude that J_i(M 2 , T, /1) —> at T — > 00. Next, from the gap equation (|22l) it follows 
that hf(T) — > at T — > 00. Therefore V(Mf,T,/i) ~ C(T, (1) at large T, i.e. we can say that the asymptotics is 
totally determined by the term C(T, /i), yielding v(T — > 00) = 31.5, which is independent of any model parameters. 
This conclusion is attractive because it agrees with the general arguments of the previous paragraph. In fact, to many 
readers our conclusion that v{T) has the correct asymptotics, may seem to be a much more compelling argument for 
fixing C(T, fx) than the assumption made that the different cut-off procedures must give the same result when cut-offs 
are removed. 

A clear insight into the origin of this result can be obtained by considering the contribution of C(T, /i) and the term 
~ { n q0 + n qo) in e q- P8"]) to the number of effective degrees of freedom f(T), because exactly these terms determine 
the correct low-temperature behavior of the function v{T). We designate this contribution by v^{T), and plot it in 
fig. [5] Thus, we consider the following function of temperature 

"A(T) = ^1°° \pe\Mpe\ (1 - PAfe)^r^ (34) 

calculated at a fixed value of the cut-off parameter A (the value taken is the same for both parameter sets considered, 
see Table I). At first glance one might wish to associate u\(T) with the contribution of massless states, introduced to 
jmed assigning the low limit M = to the integral (IM1) . However, this expectation is potentially fallacious. One 
can easily see from eq. ([34]) that v\(T) vanishes for all T if the the integral is not regularized f>Ap E ~^ !• Hence only 
the auxiliary Pauli-Villars states of mass A contribute to v\ (T) 

90N C f°° _ 4 _ ( 2 d \n qA + nq!y 
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FIG. 1: The number of effective degrees of freedom associated with Nf = 3 flavor quarks as a function of T/T c at fi = in the 
NJL model with 't Hooft and eight-quark interactions (bold solid curves correspond to our calculations with a finite cut-off, 
and thin solid curves show the same patterns when one removes the regulator in the thermal quark energy integral) compared 
with the result of Fukushima [23fl (long dashes). The critical temperature is T c — 190 MeV (finite cut-off), T c = 179 MeV 
(no regulator) for parameter set (a) (upper bold solid line, and upper thin solid line, correspondingly) , and T c = 135 MeV 
(finite cut-off), T c — 132 MeV (no regulator)for set (b) (lower bold solid line, and lower fine solid line, correspondingly). 
The long-dashed-curves are taken from fig. 3 of paper [23]: the upper one corresponds to the NJL model with 't Hooft 
interactions; the lower one to the NJL model with the Polyakov loop. The short-dashed horizontal lines correspond to the 
asymptotic high-temperature Stefan-Boltzmann (ideal gas) limit for the theory of massless fermions with three and two flavors 
respectively. 



This contribution is given by an integral with a positive integrand. Here we have two dimensionful parameters, T and 
A. Therefore, at large T the series for f\(T) can be organized in powers of the dimensionless ratio A/T, i.e. 



180 N c 



dx 



(x 2 



T 1 



1 



A 2 



2T 2 x 2 



1 



2\N r 



1 + 



(36) 



A useful observation is that although the heavy mass states determine the value of the integral, it still has the correct 
asymptotics at large T, which does not depend on A. In other words, the integral ~ (n q o + riqo) in eq. (1281) must 
vanish at T — » 00, in order that the function v(T) has the right asymptotic behavior. This really happens, due to the 
Pauli-Villars regulator pAp~ E - In fig- [Hone can see that up to around T ~ 100 MeV the contribution to the effective 
number of degrees of freedom is practically constant and nearly zero, then it increases and reaches asymptotically the 
value 31.5. The constant dashed curve shows the same quantity if the integral is not regularized. Such a big difference 
in the behavior of v\(T) translates in v(T) (see fig. [1} to an enhancement in the number of quark degrees of freedom 
around the critical temperature in the presence of a regulator. 

The conclusion that the Pauli-Villars regularization leads to the correct asymptotics is attractive because in the NJL 
model with 3-dimensional cut-off the Stefan-Boltzmann limit can be reached only when the cut-off is removed in the 
matter integrals [22j |. which are finite in themselves. Such a selective removal of the cut-off must be taken "cum grano 
salis", since technically the NJL model requires the presence of a finite ultraviolet cut-off throughout all integrals. 
Nevertheless, one might think of criticizing our result on the grounds that we seem to get the correct asymptotics due 
to the auxiliary and therefore unphysical Pauli-Villars terms. This is certainly not true. To see this let us return to 
eq. and consider now the contribution of the term ~ (n q M + w?Af) to v(T). (Our arguments will be based on 
considering the simplified case M* = = M* = M*. We need not be rigorous in the following discussion, because 
the unitary symmetry breaking effects are unimportant for the asymptotics.) The integral describes the thermal 
energy of massive quarks, but does not contribute at large T 



90N C 

7T 4 T 4 



\PE\ 4 d\p E \pAp E 



n q M, + n qM , 
E P (M*) 



2lN r 



l-l + O 



M* v^ 2 + MS 



T 



T 



(37) 



The vanishing result is a consequence of a total cancellation of two contributions: the first, 1, in the square brackets 
represents the contribution of the physical states with mass M*; the second, —1, comes from the Pauli-Villars regulator. 
What is interesting here is that the second term, if one joins it with the other unphysical contributions of eq. (|36l) . 
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FIG. 2: Solid line: degrees of freedom due to the Pauli-Villars states, v\(T), at zero chemical potential, as function of T [GeV]. 
Dashed line: the same quantity at A — Y oo. 

cancels them entirely. In other words, the correct large T asymptotics in (|33[) can equally be assigned to the pure 
physical states contribution as well. 

Let us compare our results with the recent estimates of Fukushima [23[ , made on the basis of the three- flavor N JL 
model with and without the Polyakov loop. The starting values are very similar: our curve for the set (a) agrees well 
with the Fukushima estimate made in the standard NJL model approach with the 't Hooft six-quark interactions. 
Indeed, if we remove the cut-off dependence in the quark number occupation integrals (|28[) . like it has been done 
in [23j . the curves almost coincide: compare the upper thin solid line of our set (a), where the regulator has been 
removed, with the long dashed curve of Fukushima, almost on top of it in fig. [TJ Taking systematically into account 
the finite value of the cut-off, we obtain the upper bold curve corresponding to the set (a). The set (b), compared with 
set (a), shows a rather strong (more than 50%) suppression of the abundant artificial quark excitations at T/T c > 0.2 
due to the large OZI- violating eight-quark interactions (see lower bold curve in fig. Q]for the finite cut-off result and 
lower thin solid line, where the cut-off condition has been relaxed. Notice that in both cases (a) and (b) the finite 
cut-off leads to larger values of u, because Vk{T) is non-zero and positive at A = oo; this effect is more pronounced in 
the case of set (a)). Although with set (b) the model still fails to describe accurately the pressure at 0.3 < T/T c < 1, 
it leads to an improved description of the number of quark degrees of freedom as compared to set (a). Thus the 
OZI-violating interactions could potentially play an important role in the description of quark excitations. At least 
one should not exclude the set (b) if one includes the Polyakov loop effects in the framework of effective NJL-type 
models with multi-quark interactions, see [24j for results with the 3-dimensional cut-off. 

IV. PHASE DIAGRAM 

Fig. [3] shows the phase diagrams for the two sets of parameters. One observes a larger window for the first order 
transition regime in the case of stronger 8q interaction coupling g\, the critical endpoint is situated at (he-,Te) — 
(155, 108) McV, whereas for the smaller coupling it is at (/ie,Te) — (338, 53) MeV. The temperature at zero chemical 
potential, in the crossover regime, is substantially smaller for the large 8q coupling case, around T c ~ 135 MeV 
compared to T c ~ 193 MeV for the sets (b) and (a) respectively 12]. The approximate values for T c were obtained as 
usual, through the condition c\ 2 M/AT 2 = of strongest change in the slope of the constituent quark masses. 

The length of the line of first order transitions can be increased further with increase of the 8q coupling g\ and 
after some critical value it goes all the way through to the /i = end. A corresponding set of parameters has been 
considered in [l3| . In this way we obtain the thin solid line in figs. G2 S]that connects all CEP obtained by varying 
gi in the allowed interval for stability of the effective potential. 

However we point out that also in the case of strong 8q couplings our results have with other approaches in common 
that the dynamical u and d quark masses suffer a considerable reduction at the transition temperatures, while the 
strange quark mass remains roughly the same. For instance for set (b) at /i = 170 MeV at the first order transition 
temperature T = 103 MeV, M u jumps from M u = 208 MeV to 128 MeV, while the strange quark mass changes only 
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FIG. 3: Phase diagrams indicating first order phase transition lines: the short solid line corresponds to parameter set (a) of 
weak 8q coupling constant; the long solid line is for set (b) with strong 8q coupling. Dashed lines are spinodals. Circles (blue 
online) indicate the critical endpoints (CEP) for sets (a) and (b) respectively and dotted lines correspond to the crossover 
region. The thin solid line (blue online) represent all CEP obtained by varying the 8q coupling g± and keeping as before the 
meson mass spectra at T = [i — unchanged, except for the cr-meson mass. All units are in GeV. 



Z^CEP 



FIG. 4: The line of CEP described in fig. [3] now in the ficEP,gi plane. 



from M s = 487 MeV to 448 MeV. By increasing the temperature keeping fi fixed, one obtains at T = 160 MeV the 
values M u — 24 MeV and M s — 381 MeV. Profiles of the quark masses are given in [25j for the T = Q,fi ^ case. 
The slow decrease of the strange quark mass is present in studies involving realistic values for the strange current 
quark mass [4] . 

The position of the critical endpoint in the three flavor NJL model with Ua (1) breaking has been analyzed in 
[23j |. with the parameter set of [26], ([1e,Te) = (324, 48) MeV in comparison with the case with inclusion of the 
Polyakov loop, (/xb,Te) = (313, 102) MeV. The increase of Te is explained by the suppression of the artificial quark 
excitations at finite temperature and density in the presence of the Polyakov loop contribution f23j j . Our set (a) with 
(he,Te) — (338, 53) MeV is in reasonable agreement with the above estimates. 
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Our result for set (b), {he,Te) — (155, 108) MeV, yields a twice lower value for the critical he as compared to 
set (a). This new feature is related with the OZI- violating eight-quark interactions: it is well-known that NJL-type 
models without 8j£-forces have the tendency in common to lead to a critical endpoint at relatively high density, above 
[j,E ~ 300 MeV [23[ . The small value found here is a good indicator that the main force responsible for dynamical 
chiral symmetry breaking has changed, being now associated with the 't Hooft 6q-interactions. 

Let us also take notice of the difference between the critical temperatures for the case (b), T c = 108 MeV, and the 
Fukushima's case with Polyakov loop, T c = 204.8 MeV. We expect that the inclusion of the Polyakov loop in our 
analysis will increase somewhat the value of T c . 

The NJL model with 8q interactions and without UaO-) symmetry breaking has been analyzed as well for the SU(2) 
flavor case [13, l28j . Although the 8q forces are not needed to stabilize the effective potential of the model in the 
two flavor case, the same tendency as above in i) was observed, for instance: (he,T e ) = (276, 62) MeV with the 8q 
interactions, whereas (/i E ,T E ) — (330, 47) MeV without them (27j |. 

Another variant without 8q terms of the two flavor NJL model, includin g th e vector-isoscalar interactions, which 
induce an effective chemical potential, has been considered a long time ago [29| . In this case, the critical endpoint is 
located at larger \i and lower temperatures (fi E ,T E ) = (350, 40) MeV. This effect has been studied also recently in 

MM- 

Further results concerning the position of the critical endpoint within other model approaches is given in the review 
Q (note the notation there is in terms of hb = 3/i), in comparison with lattice results and the freezeout points 
extracted from heavy-ion experiments. For our two sets the critical endpoints are situated in case (a) slightly above 
and in case (b) slightly below the freezeout points obtained at different collision energies [H, [3(| ■ 

V. CONCLUSIONS 

The thermodynamic potential of the three flavor NJL model with 't Hooft and 8q interactions has been obtained 
in stationary phase approximation (at leading order) and using the Pauli-Villars regularization in quark loops (at 
one-loop level). 

The main conclusions of this work can be classified as follows. 

Firstly, we argued that a non-renormalizable effective theory at finite temperature and chemical potential, can 
be self-consistently studied with the use of the Pauli-Villars regularization. The importance of the regulator in the 
matter parts for consistency has also been discussed in [32j in connection with correlators. Instead, our present study 
is devoted to the construction of the thermodinamic potential. Unlike the conventional approach with 3-dimensional 
cut-off, the Pauli-Villars technique leads to the right asymptotic behavior of relevant thermodynamic observables. It 
is one of the main results of this paper. 

Secondly, we quantified the effect of the new 8q terms on the number of degrees of freedom and on the phase 
diagram. We observe that in the large 8q coupling regime a strong depletion of the number of degrees of freedom 
below T c is reached in comparison with the weak coupling case, working in the same direction as the effect produced 
through inclusion of the Polyakov loop in the model without 8q interactions. 

Thirdly, we conclude that the NJL model with the 8q stabilizing interactions does not impede the possibility of 
having a phase diagram consisting only of first order transitions even for realistic quark masses. This will depend on 
the strength of the OZI-violating 8q interactions. At /i = there is growing evidence from lattice calculations that the 
transition is a crossover [3.1] . This would set an upper limit for the 8q coupling, which nevertheless can be sufficiently 
strong to trigger the first order transitions regime at low values of h ^= 0. This point deserves to be studied more 
carefully especially because it can help us to clarify the quark dynamics which is responsible for the mechanism of 
spontaneous chiral symmetry breaking. Indeed as it is shown above the shift of the critical endpoint to lower values 
of he is possible only when the 6q interactions are responsible for the chiral phase transition (set (b)); however if the 
4<7 coupling G exceeds its critical value, the 4q interactions drive the chiral phase transition in NJL-type models and 
as a result the critical endpoint is located at large values of he (set (a)). 
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